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AN  EMPIRICAL  MODEL  TO  COMPUTE  THE 
VELOCITY  HISTORIES  OF  FLYERS 
DRIVEN  BY  ELECTRICALLY  EXPLODING  FOILS 


by 

S.  C.  Sriimidt,  W.  L.  Seitz,  and  Jerry  Wackerle 
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A  modification  of  the  gas  gun  (Gurney)  formulation  is  used  to  compute  the 
velocity  and  position  histories  of  flyers  driven  by  electrically  exploded  metal 
foils.  The  model  is  based  on  a  numerics.1  time  integration  of  an  energy  con¬ 
servation  statement  for  the  flyer  and  the  expanding  high-pressure  metal 
vapor.  Empirically  altered,  experimental  power  curves  are  used  for  the 
time-dependent  energy  term  in  the  conservation  equation.  Computed  burst 
times  and  flyer  velocity  histories  for  1.5-  to  25-mm-square  aluminum  foils 
agree  favorably  with  available  experim  ental  data.  Comparison  of  calculated 
and  measured  results  for  single  cases  with  exploding  copper  and  magnesium 
foils  suggests  that  the  model  is  also  applicable  for  these  materials. 


I.  INTRODUCTION 

During  the  past  decade,  electrically  driven  flyer 
systems  have  evolved  into  effective  laboratory  high- 
pressure  shock-wave  generators,  which  are  par¬ 
ticularly  useful  for  investigations  of  the  short- 
duration  shock  initiation  of  condensed  explosives.'  A 
typical  arrangement  for  accelerating  thin  flyers  is 
shown  in  Fig.  1.  A  capacitor  bank  is  discharged 
through  a  thin  metal  foil,  effecting  its  abrupt 
vaporization.  This  "burst"  of  the  foil  generally  occurs 
with  a  sharp  maximum  foil  resistance  and  electrical 
power  input.  The  plastic  layer  is  sheared  at  the  inner 
radius  of  the  barrel,  and  the  flyer  disk  so  formed  is 
subsequently  accelerated  by  the  high-pressure 
vapor.  Flyer  velocities  as  high  as  14  mm/MS  have 
been  achieved  in  this  way.' 

The  gas  gun  (Gurney)  approximation  is  a  simple 
analysis  that  has  been  used  to  predict  the  motion  of 
plates  and  shells  driven  by  detonated  chemical 
explosives,*^  and  recently  has  been  used  to  treat 
electrically  driven  flyers.'-*""  Conceptually,  the 
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physical  system  modeled  is  the  expansion  of  high- 
pressure  gas,  initially  at  rest,  against  confining 
plates  or  shells.  A  "Gurney  energy"  is  assumed  to  be 
deposited  uniformly  throughout  the  gas  before  any 
motion  begins.  This  energy  is  taken  as  an  empirical 
fraction  of  the  energy  of  reaction  for  chemical 
explosives*'*  and  has  been  correlated  to  the  joule 
heating  near  burst  for  electrically  exploded  foils. 

For  explosively  driven  flyers,  the  gas-gun  ap¬ 
proach  predicts  flyer-velocity  histories  that  agree 
favorably  with  both  experimental  results  and 
hydrodynamic  computer  code  calculations.*’* 
Previous  application*’*'"  of  the  Gurney  formulation 
to  electrically  driven  flyers  has  been  limited  to  the 
calculation  of  terminal  velocities,  and  has  not 
produced  entirely  satisfactory  results.  The  standard 
gas-gun  formulation  does  indeed  predict  that  ter¬ 
minal  velocity  is  achieved  in  typical  run  distances; 
however,  a  much  more  gradual  and  continuing  flyer 
acceleration  is  observed.  Since  the  corresponding  ex¬ 
ploding  foil  power  histories  measured  by  us  and 
others*’"  suggest  additional  energy  deposition  in  the 
metallic  vapor  subsequent  to  foil  burst,  a  formula¬ 
tion  of  the  treatment  with  time-dependent  energy 
deposition  is  indicated. 

In  the  work  described  here,  a  conceptually 
satisfactory  mod:  flcation  of  the  gas-gun  treatment  is 
used  to  compute  velocity  histories  for  flyers  driven 
by  electrically  exploded  foils.  Rather  than  assume 
instantaneous  energy  deposition,  we  empirically 
relate  the  joule  heating  of  the  foil  material  to  the 
observed  electrical  power  histories,  with  allowance 
for  energy  dissipated  in  forming  the  metal  plasma. 
These  changes  yield  computed  velocity  histories 
considerably  different  from  those  resulting  from  the 
usual  Gurney  model  and  in  favorable  agreement 
with  available  experimental  data.  In  many  in¬ 
stances,  the  model  also  affords  reasonable  estimates 
of  foil  burst  time. 

The  model  is  strictly  empirical.  Even  though 
many  of  the  relevant  physical  processes  that  occur 
during  foil  heating  and  expansion  are  used  to  incor¬ 
porate  the  empiricisms,  many  more  processes  (e.g., 
plasma  recombination,  spatially  nonuniform  energy 
deposition)  are  ignored.  The  principal  improvement 
over  previous  models  of  electrically  driven  flyers  is 
that  a  computation  of  complete  velocity  histories 
can  be  made  for  a  wide  range  of  foil  and  flyer  dimen¬ 
sions. 


II.  ANALYSIS 


A.  Simple  Gas-Gun  Formulation 


The  usual  Gurney  model  is  formulated  by  the  con¬ 
servation  of  the  total  energy  (internal  and  kinetic)  of 
the  driving  gas-flyer  system.  Initially,  the  system  is 
at  rest  and,  with  compression  and  internal  energy  of 
the  flyer  neglected,  the  total  energy  is  the  internal 
energy  of  the  gas.  This  Gurney  energy.  Eg,  is  as¬ 
sumed  to  be  generated  instantaneously  from  an  ex¬ 
plosive  reaction  or  from  electrical  heating  of  a  foil 
vapor.  In  the  subsequent  expansion  of  the  gas,  the 
internal  energy,  E,  density,  p,  and  pressure,  p,  are 
taken  to  be  spatially  uniform.  The  mass  velocity,  u, 
is  assumed  to  vary  linearly  with  distance,  r,  in  the 
gas  and  to  have  the  uniform  value  ur  in  the  flyer. 

For  electrically  driven  flyers,  slab  geometry  with 
the  vapor  sandwiched  between  the  flyer  and  an  in¬ 
finitely  massive  tamper  is  assumed.  With  the  above 
assumptions,  the  energy  conservation  statement  is 
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where  mg  and  mr  are  the  gas  and  flyer  masses  per 
unit  area  and  rt  the  position  of  the  gas-flyer  inter¬ 
face.  With  u(r,t)  =  (r/rf(t))Uf(t)  (imd  noting  that 
p(t)rf(t)  =  mg),  this  expression  reduces  to 


(2) 


where  R  =  mr/mg.  Assuming  that  the  equation  of 
state  of  the  vapor  is  that  of  a  polytropic  gas. 
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where  the  gas  constant,  y  =  5/3,  is  chosen  on  the  as¬ 
sumption  that  the  metal  is  monatomic.  When  the 
pressure  is  considered  as  an  accelerating  force  ex¬ 
erted  per  unit  area  of  the  flyer, 


p 


(4) 
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the  internal  energy  can  be  eliminated  from  Eq.  (2), 
yielding 


(5) 


The  transformation  duf/dt  =  l/2(dUfVdrf)  converts 
Eq.  (5)  to  a  first-order  differential  equation  for 
Uf’{rf),  which  can  be  integrated  from  the  initial  gas- 
flyer  interface  position,  r*,  to  an  arbitrary  final  posi¬ 
tion,  rr,  to  give 
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where  ^  =  2/3  2/(9<S.  As  the  flight  distance 

becomes  much  larger  than  the  foil  thickness,  rt  »  rg, 
a  terminal  velocity 
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the  simple  analysis  is  applied  to  the  other  flyer/foil 
assemblies  described  here. 

The  failure  of  the  flyer  to  reach  a  terminal  velocity 
as  quickly  as  predicted  can  be  understood  by  ex¬ 
amining  the  power  curves  for  the  exploding  foils, 
such  as  curve  A  in  Fig.  2.  Contrary  to  the  instan¬ 
taneous  energy  deposition  at  burst,  which  is  as¬ 
sumed  when  using  the  Gurney  model,  the  actual 
power  deposited  is  not  sharply  peaked  near  the  time 
of  burst,  but  exhibits  a  moderately  fast  rise  to  peak 
power  and  a  more  gradual  decay  after  burst.  The 
gradual  decay  suggests  that  joule  heating  of  the  foil 
is  still  occurring  at  relatively  large  postburst  times. 
Addition  of  energy  to  the  expanding  gas  would  sus¬ 
tain  a  greater  gas  pressure  than  expected  from  the 
simple  model,  and  would  lead  to  a  continuing  ac¬ 
celeration  of  the  flyer.  Thus  the  energy  conservation 
statement  should  be  formulated  with  a  time- 
dependent  energy  deposition  term  to  allow  closer 
agreement  with  experimental  powei  curves. 


is  attained. 

As  mentioned  earlier,  previous  applications  of  the 
Gurney  method  to  electrically  driven  flyers  have  all 
been  made  on  the  basis  of  a  terminal  velocity  such 
as  that  expressed  in  Eq.  (7 ).*•*■"  The  Gurney  energy 
was  correlated  to  the  burst  current  density,  Jg,  in  the 
form 

Eg  -  KJb  *  (8) 

where  K  and  n  were  empirically  determined  cons¬ 
tants  for  a  gi>'en  foil  material  and  for  velocity 
measurements  at  a  specified  flight  distance.  This 
model  successfully  predicts  the  flyer  velocity 
dependence  on  the  burst  current  density  and  flyer 
thickness  and  is  adequate  if  these  are  the  only 
parameters  to  be  varied. 

The  consideration  of  only  the  terminal  velocity  for 
electrically  driA'en  flyers  is  consistent  with  the  as¬ 
sumption  of  the  simple  model  For  typical  con¬ 
figurations  and  flight  distances  of  interest,  the  input 
parameters  to  Eq.  (6)  would  be  0  s  1  and  rg/ri  3 
0.01,  and  terminal  velocity  is  nearly  achieved. 
However,  under  these  conditions  terminal  velocities 
are  not  usually  observed.  For  example,  the  velocity 
history  measurement  of  Weingart  and  coworkers* 
(Fig.  2,  curve  B,  of  the  present  report)  has  a  much 
more  gradual  acceleration  than  the  trajectory,  B', 
calculated  with  the  Gurney  method  and  the  burst- 
current  correlation.  Similar  comparisons  occur  when 


B.  Time-Dependent  Energy  Deposition 


If  the  constant  Gurney  energy.  Eg,  is  replaced  by  a 
time-dependent  deposited  specific  energy,  Q(t),  the 
conservation  statement,  Eq.  (5),  can  be  rewritten  as 


^  [-(3  *•)  W  *  »">]  ■ 

Initial  conditions  for  the  integration  are  Q(0)  =  ur(0) 
=  0  and  rf(0)  =  r,.  With  Q(t)  specified  by  observed 
power  histories,  a  numerical  solution  is  necessary. 
Specifically,  the  set  of  three  coupled  differential 
equations,  Eq.  (9). 


drj(t) 

dc 


-  Uj(t) 


(10) 


and 


P(t)  (11) 

are  integrated  numerically  using  LASL  TLIB 
Subroutine  ODE.  The  power  history,  P(t),  is  derived 
from  a  separate  analysis  of  foil  current  and  voltage 
measurements  and  a  fifth-order  polynominal  inter¬ 
polation  scheme  (LASL  TUB  Subroutine  AK- 
NINT).  Initial  calculations  indicated  a  need  to 
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Time  (^s) 

Fig.  2. 

Flyer  velocity  and  foil  power  va  time;  0.061 -mm- thick,  2S.4~mm-8quare  aluminum  foil;  0.25- 
mm-thick,  25.4-mm-diam  Mylar  flyer;  40-kV  discharge  voltage.  Curve  A,  experimental 
power  history.  Curve  B.  experimental  velocity  history.  Curve  B',  computed  velocity  history 
using  Gurney  model.  Curve  C,  computed  velocity  history  using  unaltered  postburst  power 
history.  Curve  D,  computed  velocity  history  using  modified  postburst  power  history. 


reduce  the  actual  power  input.  The  reduction  is  ac¬ 
complished  by  empirical  fitting. 

III.  EMPIRICAL  FITTING 

The  analysis  is  carried  out  with  the  aid  of  two  em¬ 
piricisms:  (1)  a  preburst  energy,  I,  empirically 
chosen  as  the  sum  of  the  heats  of  fusion,  vaporiza¬ 
tion,  and  ionization,  is  subtracted  from  the  ex¬ 
perimental  power  curve  so  that  Q(t)  -  0  in  Eq.  (8) 
until  the  /  P(t)dt  =  I,  and  (2)  modification  of  the 
poatburst  portion  of  the  power  curve  so  that 
predicted  velocity  histories  agree  with  measure¬ 
ments. 


A.  Energy  Used  for  Fusion,  Vaporization,  and 
Ionization 

The  first  empiricism  is  suggested  from  experimen¬ 
tal  velocity  and  power  histories.  As  noted  in  Fig.  2, 
flyer  movement  starts  very  close  to  the  time  of  peak 
power  or  foil  burst,  implying  that  energy  deposited 
in  the  foil  before  burst  does  not  contribute  directly 
to  flyer  acceleration.  Foil  melting,  vaporization,  and 
ionization  are  possible  dissipative  processes  to  ab¬ 
sorb  this  energy.*  As  a  first  approximation,  I  was  set 
equal  to  the  sum  of  the  heats  of  fusion,  vaporization. 


*The  actual  heating  path  may  not  include  a  preaaure-volume 
atate  which  allows  for  a  vaporization  transition. 
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and  ionization  at  atmospheric  pressure;  values  for 
the  three  foil  materials  of  interest  are  given  in  Table 
I. 

Calculations  were  performed  for  experiments  by 
Weingart  and  coworkers*  and  by  Stanton**  by  using 
their  measured  power  curves  and  the  I  =  32.3  J/mg 
value  for  aluminum,  with  the  results  given  in  Figs.  2 
and  3.  In  both  cases,  32.3  J/mg  corresponded  to  the 
integration  of  the  specific  power  (curve  A)  to  a  time 
near  burst.  In  both  instances,  the  calculated  velocity 
histories  (labeled  C)  gave  times  of  flyer  movement 
coinciding  with  both  the  times  of  initial  motion 
observed  (curves  B)  and  the  times  of  peak  power  in 
the  foils.  Following  burst,  the  predicted  flyer 
velocity  histories  have  the  correct  shape,  but  are 
considerably  larger  than  those  observed.  This  dis¬ 
agreement  dictates  a  further  modification  of  power 
input,  P(t). 

B.  Modified  Power  History 

Even  though  the  correction  is  empirical,  some 
rationale  for  a  modification  of  the  measured  power 
histories  can  be  argued.  Following  burst,  the 
material  accelerating  the  flyer  is  a  highly  ionized, 
highly  conductive  plasma;  however,  a  region  of 
higher  resistance  liquid  and  un-ionized  vapor  must 
exist  between  the  plasma  and  the  solid  lead  (see  Fig. 
4).  Because  voltage  probes  must  be  positioned  so 
that  the  melt/vapor  transition  region  is  included  in 
the  measurement,  the  electrical  resistance  and 
power  observed  do  not  properly  reflect  the  energy 
dissipated  in  the  plasma. 


Another  source  of  disagreement  between  theory 
and  experiment  is  the  assumption  of  an  infinitely 
massive  tamper  in  the  development  of  Eq.  (9)  so 
that  the  kinetic  energy  (and  velocity)  of  the  flyer  is 
not  diminished  by  the  kinetic  energy  of  tamper 
material.  This  assumption  is  unrealistic,  particular¬ 
ly  during  the  period  just  following  burst,  when  pres¬ 
sures  of  several  gigapascals  are  estimated  in  the  foil 
plasma.  A  calculation  to  correct  for  this  effect  was 
made  with  a  "growing  tamper"  model,  in  which  the 
mass  per  unit  area  of  the  confining  tamper  material 
was  assumed  to  be  zero  at  burst  and  to  subsequently 
increase  according  to  its  density  multiplied  by  its 
characteristic  sound  speed.  Appropriate  modifica¬ 
tion  of  the  energy  conservation  relation  and  explicit 
use  of  momentum  conservation  gave  a  tractable 
problem,  and  calculated  velocity  histories  had  the 
expected  improved  agreement  with  observation. 
Still,  the  need  for  a  large  empirical  correction 
remained;  consequently,  we  included  the  effect  of 
tamper  motion  in  a  correction  rather  than  employ¬ 
ing  the  more  complicated  formulation  of  the  growing 
tamper  model. 

A  correction  factor  applied  to  the  postburst  power 
density  was  chosen  to  force  agreeme  nt  for  cases  with 
25.4-  and  9.53-mm-8quare  aluminum  foils  (Figs.  2 
and  3).  The  form  used  was 

Correction  Factor  »  [0.958  -  0.166  U) 


Here  W  is  the  width  of  the  foil  in  millimeters,  T  is 
the  time  from  foil  burst  (calculated  as  described 
above),  and  6  is  chosen  as  0.1  ms  for  the  large  foils 


TABLE  I 

MATERIAL  PROPERTIES  OF  FOIL  METALS 
(Taken  from  Ref.  12) 


Material 

Density 

(mg/mm*) 

Heat  of 
Fusion 
(J/mg) 

Heat  of 
Vaporization 
(J/mg) 

Heat  of 
Ionization 
(J/mg) 

Aluminum 

2.79 

0.40 

10.5 

21.4 

Magnesium 

1.74 

0.37 

5.42 

30.3 

Copper 

8.89 

0.21 

4.80 

11.7 
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Flyer  velocity  and  foil  power  vs  time;  0.051 -mm-thick,  9.53-mm-square  aluminum  foil;  0. 127- 
mm-thick,  9.53-mm-diam  Mylar  flyer;  18-kV  discharge  voltage.  Curve  A,  experimental 
power  history.  Curve  B,  experimental  velocity  history.  Curve  C,  computed  velocity  history 
using  unaltered  postburst  power  history.  Cur, re  D,  computed  velocity  history  using  modified 
postburst  power  history. 


Melt  and  vapor  region  of  exploding  foil. 


under  discussion  and  0.025  iis  for  the  smaller  foils 
described  later.  The  results  of  applying  this  post- 
burst  correction  factor  are  shown  by  curves  D  in 
Figs.  2  and  3.  Although  the  results  are  superior  to 
those  obtained  with  the  simple  Gurney  formulation, 
the  real  test  of  the  model  is  its  ability  to  simulate 
velocity  histories  for  foil/flyer  geometries  much  dif¬ 
ferent  from  those  used  to  calibrate  the  empirical  cor¬ 
rection  factor. 


IV.  COMPARISON  OF  MODEL  WITH  AD¬ 
DITIONAL  EXPERIMENTAL  DATA 

In  this  section  the  empirical  fit  structured  above  is 
applied  to  a  variety  of  different  foil/flyer  configura¬ 
tions  to  validate  its  usefulness  as  a  predictive  tool 
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and  also  to  determine,  as  well  as  possible,  its  range 
of  applicability. 

A.  Large  Aluminum  Foil  Data 

Weingart  and  coworkers*  **  have  studied  flyers 
driven  by  aluminum  foils  of  different  dimensions 
and  exploded  with  different  firing-set  voltages. 
Figures  5A  and  5B  give  an  experimental  power 
history  and  a  single  velocity  datum**  for  a  foil/flyer 
configuration  and  discharge  voltage  similar  to  that 
given  in  Fig.  2,  except  that  the  aluminum  foil 
thickness  has  been  increased  by  a  factor  of  3.  Also 
shown  are  the  predicted  velocity  history  and  the 
velocity-distance  dependence  computed  using  the 
empirical  fit  discussed  above.  Both  the  computed 
velocity  and  the  estimated  burst  time  agree 
favorably  with  experimental  results.  Even  though 
the  exploding  foil  for  this  case  has  tripled  in 
thickness  from  the  example  shown  in  Fig.  2,  the 
postburst  velocity  histories  of  the  two  configurations 
are  comparable.  The  specific  power  in  the  thinner 
foil  is  considerably  greater  than  in  the  thicker  foil, 
but  the  total  energy  deposited  is  approximately  the 
same  for  the  two  configurations.  The  thicker  foil  re¬ 
quires  more  time  and  three  times  the  electrical 
energy  to  burst,  but  since  1  is  a  relatively  small  frac¬ 
tion  of  the  total  electrical  input  energy,  the 
velocities  attained  are  comparable.  This  again 
emphasizes  the  importance  of  considering  the  post¬ 
burst  power  input. 

Figures  6A  and  6B  show  a  power  history  and  single 
velocity**  datum  for  a  foil/flyer  configuration  in 
which  the  foil  width  is  half  that  of  the  configuration 
of  Fig.  2  and  the  capacitor  discharge  voltage  is  also 
decreased  by  the  same  factor.  As  can  be  seen  from 
the  predicted  velocity  history,  a  good  estimate  of  the 
foil  burst  time  is  obtained  because  initial  flyer 
movement  occurs  at  approximately  the  time  of  peak 
power.  Figure  6B  also  shows  that  the  velocity  com¬ 
puted  using  the  empirical  model  agrees  favorably 
with  the  measured  velocity. 

B.  Large  Magnesium  and  Copper  Foils 

To  assess  the  model's  range  of  validity,  it  was 
compared  with  Stanton's**  observations  of  flyers 


driven  by  electrically  exploded  magnesium  and  cop¬ 
per  foils  with  geometries  and  firing-set  conditions 
similar  to  those  of  the  aluminum  foil  experiment 
(Fig.  3).  The  empirical  postburst  specific  power  cor¬ 
rection,  calibrated  using  the  aluminum  foil  observa¬ 
tions,  was  applied  to  the  cases  of  magnesium  and 
copper  foils. 

The  predicted  velocity  history  for  the  magnesium 
foil  (Fig.  7)  is  in  reasonable  agreement  with  obser¬ 
vation — the  small  discrepancy  being  due  mainly  to 
the  calculated  later  start  for  the  flyer  motion.  This 
disagreement  occurs  despite  the  fact  that  the  com¬ 
puted  burst  time  corresponds  to  the  peak  in  the 
power  history. 

For  the  copper  foil,  the  computed  flyer  velocity 
history  (Fig.  8,  curve  C),  predicts  a  burst  time  and 
initial  flyer  motion  that  occur  considerably  later 
than  the  peak  in  the  power  history  or  initial 
measured  flyer  movement.  If  the  ionization  energy 
of  copper  is  not  included  in  the  determination  of  I, 
the  computed  velocity  history  (curve  D)  agrees  with 
that  observed.  Conceivably,  the  lower  specific  power 
might  provide  some  argument  for  assuming  in¬ 
complete  ionization  of  the  foil;  however,  this  deduc¬ 
tion  should  be  accepted  with  caution  because  of  the 
highly  empirical  character  of  the  analysis. 

C.  Small  Aluminum  Foil  Data 

We  have  measured  velocity,  current,  and  voltage 
histories  for  seven  1.52-mm-diam  Mylar  flyers 
driven  by  0.011-mm-thick,  1.52-mm-square 
aluminum  foils.  A  2-^F  capacitor  discharge  unit, 
charged  to  3  kV,was  used  to  explode  the  foils.  Four 
tests  were  fired  with  0.75-mm-long,  1.52-mm-diam 
barrels  attached  to  the  Mylar  surface  (see  Fig.  1). 
Velocities  at  the  end  of  the  barrels  were  determined 
from  streak-camera  records  of  the  flyers  striking 
~0.15-mm  Lucite  step  flashers.  The  other  three  ex¬ 
periments  were  performed  without  barrels,  with 
velocity  histories  obtained  using  a  modification  of 
the  streak-camera  reflection  technique.**  Flyer  mo¬ 
tion  was  recorded  by  observing,  at  an  oblique  angle, 
the  light  emitted  from  the  exploding  foil  and  trans¬ 
mitted  through  uncoated  slits  on  the  flyer  surface. 
Barrels  were  not  used  in  these  experiments  because 
they  would  partially  obscure  the  flyer  surface  at 
small  travel  distances.  In  the  absence  of  a  fiducial 
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Fig.  5. 

Flyer  velocity  and  foil  power  va  time  (Fig.  5A)  and  flyer  velocity  va  travel  diatance  (Fig.  5B); 
0.l52-mm-thick,  25.4-mm-aquare  aluminum  foila;  0.25-mm-thick,  25.4-mm-diam  Mylar 
flyera;  40~kV  discharge  voltage. 


Distance  (mm) 


Fig.  6. 

Flyer  velocity  and  foil  power  V8  time  (Fig.  6A)  and  flyer  velocity  vs  travel  distance  (Fig.  6B); 
0.051 -mm-thick,  12.7-mm-square  aluminum  foils;  0.25-mm-thick,  12.7-mm-diam  Mylar 
flyers;  20.4-kV  discharge  voltage. 
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Fig.  7. 

Flyer  velocity  and  foil  power  vs  time;  0.051 -mm-thick,  6.35-mm-square  magnesium  foil; 
0.127-mm-thick,  6.35-mm-diam  Mylar  flyer;  18-kV  discharge  voltage. 
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Fig.  8. 

Flyer  velocity  and  foil  power  vs  time;  0.036-mm-thick,  9.53-mm-square  copper  foil;  0.127- 
mm-thick,  9.53-mm-diam  Mylar  flyer;  18-kV  discharge  voltage.  Curve  A,  experimental 
power  history.  Curve  B,  experimental  velocity  history.  Curve  C,  computed  velocity  history 
including  ionization  energy.  Curve  D,  computed  velocity  history  not  including  ionization 
energy. 
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on  the  streak-camera  film,  the  electrical  current 
and  voltage  records  were  related  to  the  optical  ob¬ 
servations  by  correlating  the  maximum  in  the 
voltage  measurement  with  the  first  observation  of 
light  from  the  bursting  foil. 

The  velocity  and  specific  power  histories  and  the 
change  of  velocity  with  travel  distance  deduced  from 
the  above  data  for  0.025-,  0.051-,  and  0.127-mm- 
thick  flyers  are  given  in  Figs.  9  and  10.  The  different 
velocity  curves  correspond  to  the  measurements  at 
different  points  on  the  flyer  surface.  The  initial 
velocity  of  each  of  the  flyers  corresponds  to  the 
measured  velocity  at  the  time  su^icient  light  is  first 
observed  from  the  bursting  foil.  Because  all  three 
unconflned  flyers  were  Hrst  observed  with  a  nonzero 
initial  velocity,  these  histories  have  been  shifted 
0.05  MS  to  the  right  of  estimated  foil  burst  time.  The 
estimated  error  for  the  velocity  measurements  is  ap¬ 
proximately  ±10%,  except  for  small  travel  times 
and  distances  which  have  somewhat  larger  errors. 

The  calculated  velocity  histories  and  velocity- 
distance  relationships  agree  with  observations  to 
within  experimental  error,  even  though  the  calibra¬ 
tion  of  the  specific  power  correction  factor  was  made 
using  much  larger  foil  configurations.  The  agree¬ 
ment  for  the  small  systems  could  be  improved 
further  by  altering  the  empirical  correction  factor  to 
match  the  small-foil  data. 

The  step-flasher  measurements  taken  with  barrel 
assemblies  agree  within  experimental  error  with  the 
reflection  technique  observations  taken  without  bar¬ 
rels  and  with  the  computed  velocities.  However,  the 
suggestion  that  the  confining  effect  of  barrels  is  un¬ 
important  should  be  accepted  with  caution. 

Because  the  specific  power  histories  in  Fig.  9  differ 
slightly,  the  flyer  thickness,  d,  is  the  principal 
parameter  varied  in  this  series  of  experiments. 
Observed  late-time  velocities  are  roughly  propor¬ 
tional  to  vTTd,  so  that  the  simple  gas-gun  model,  us¬ 


ing  a  correctly  selected  Gurney  energy,  would  ap¬ 
pear  to  provide  a  decent  approximation  for  this  par¬ 
ticular  parameter  variation.  However,  terminal 
velocities  computed  with  Eqs.  (7)  and  (8),  using 
measured  burst  currents  and  constants  calibrated  to 
large  aluminum  foils,  are  larger  by  a  factor  of  2  than 
those  observed  experimentally. 


V.  SUMMARY 

A  version  of  the  gas-gun  (Gurney)  model  that  al¬ 
lows  continuous  electrical  energy  deposition  has 
been  developed  and  used  to  calculate  velocity 
histories  of  flyers  driven  by  electrically  exploded 
foils.  The  time-dependent  energy  input  is  related  to 
observed  power  histories,  with  a  semi-empirical 
treatment  of  the  energy  required  to  form  the  foil 
material  plasma  and  a  purely  empirical  correction  of 
the  postburst  contribution.  Agreement  with  experi¬ 
ment  is  obtained  for  aluminum-foil  systems  varying 
by  an  order  of  magnitude  in  foil  and  flyer  dimen¬ 
sions  and  for  assemblies  with  magnesium  and  cop¬ 
per  foils.  Presumably,  the  tested  range  of  ap¬ 
plicability  of  the  model  can  be  extended  to  include  a 
greater  variety  of  geometrical  configurations  and 
other  flyer  and  foil  materials. 
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Fig.  9. 

Flyer  velocity  and  foil  power  vs  time;  0.011-mm-thick,  1.524-mm-square  aluminum  foils; 
1.524-mm-diam  Mylar  flyers;  ~3-kV  discharge  voltage. 

Figure  9A,0.025-mm-thick  flyer;  open  circle  indicates  step  flasher  data. 

Figure  9B,  0.051-mm-thick  flyer;  open  squares  indicate  step  flasher  data. 

Figure  9C,  0.127-mm-thick  flyer;  open  triangle  indicates  step  flasher  data. 
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Fig.  10. 

Flyer  velocity  vs  travel  distance;  O.Oll-mm-thick,  1.524-mm-square  aluminum  foils;  1.524- 
mm-diam  Mylar  flyers;  3-kV  discharge  voltage.  Curves  A,  0.025-mm- thick  flyer;  solid  circle, 
computed  velocity;  open  circle,  step  flasher  data.  Curves  B,  0.051 -mm-thick  flyer;  solid 
squares,  computed  velocity;  open  squares,  step  flasher  data.  Curves  C,  0.127-mm-thick  flyer; 
solid  triangles,  computed  velocity;  open  triartgle,  step  flasher  data. 


REFERENCES 

1.  R.  C.  Weingart,  R.  S.  Lee,  R.  K.  Jackson,  and 
N.  L.  Parker,  "Acceleration  of  Thin  Flyers  by  Ex¬ 
ploding  Metal  Foils;  Application  to  Initiation 
Studies,"  Sixth  Symposium  (International)  on 
Detonation,  Preprints,  p.  201  (August  1976). 

2.  R.  W.  Gurney,  "The  Initial  Velocity  of  Frag¬ 
ments  from  Bombs,  Shells  and  Grenades,”  Ballistic 
Research  Laboratory  report  No.  405  (September 
1943). 

3.  I.  G.  Henry,  "The  Gurney  Formula  and  Related 
Approximations  for  High  Explosive  Deployment  of 
Fragments,"  Hughes  Aircraft  Corp.  report  PUB-189 
(April  1967). 


4.  M.  Defoumeaux,  'Transport  of  Energy  in  Com¬ 
bustions  and  Detonations  with  Confinement,” 
Astronaut.  Acta  17,  609  (1972). 

5.  J.  E.  Kennedy  and  A.  C.  Schwijz,  "Detonation 
Transfer  by  Flyer  Plate  Impact,"  Eighth  Sym¬ 
posium  on  Explosives  and  Pyrotechnics,  The 
Franklin  Institute  Research  Laboratories, 
Philadelphia,  PA  (February  1974). 

6.  J.  E.  Kennedy,  "Gurney  Energy  of  Explosives: 
Estimation  of  the  Velocity  and  Impulse  Imparted  to 
Driven  Metal,"  Sandia  Laboratories  report  SC-RR- 
70-790  (1970). 


13 


7.  S.  J.  Jacobs,  "The  Gumey  Formula:  Variations 
on  a  Theme  by  Lagrange,"  Naval  Ordnance 
Laboratory  report  No.  74-86  (June  1974). 

8.  J.  Wackerle,  "Applied  Initiation  Studies:  Sim¬ 
ple  Considerations  for  the  Design  of  Flying-Plate 
Initiation  Systems,"  LASL  Group  WX-7  internal 
report  (October  1974). 

9.  T.  J.  Tucker  and  P.  L.  Stanton,  "Electrical 
Gumey  Energy:  A  New  Concept  in  Modeling  of 
Energy  Transfer  from  Electrically  Exploded  Con¬ 
ductors,"  Sandia  Laboratories  report  SAND-75-0244 
(1975). 

10.  T.  J.  Tucker  and  R.  P.  Toth,  'The  Electrical 
Behavior  of  Small  Exploding  Foil  Flyer  Assemblies," 
Sandia  Laboratories  report  SAND-75-0445  (1975). 


11.  P.  L.  Stanton,  The  Acceleration  of  Flyer  Plates 
by  Electrically  Exploded  Foils,"  Sandia 
Laboratories  report  SAND-75-0221  (1976). 

12.  Handbook  of  Chemistry  and  Physics,  51st  Ed. 
(TTie  Chemical  Rubber  Company,  Cleveland,  OH, 
1970). 

13.  R.  K.  Jackson,  Lawrence  Livermore  Laboratory, 
personal  communication  (November  1976). 

14.  W.  C.  Davis  and  B.  G.  Craig,  "Smear  Camera 
Technique  for  Free-Surface  Velocity  Measurement," 
Rev.  Sci.  Instrum.  32,  579  (1961). 


14 


